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Abstract 


Two  computer-based  methodologies  which  make  use  of  the  spatial  informa¬ 
tion  contained  in  the  Conservation  Needs  Inventory  (CNI)  are  presented  and 
illustrated  with  applications  to  the  Chowan-Pasquotank  River  Basin.  The 
amount  and  location  of  land  clearing  and  the  amount  and  location  of  gross 
soil  erosion  from  cropland  are  determined  and  displayed  on  computer-generated 
maps.  Potential  appl i cations  to  river  basin  planning  and  other  studies  are 
discussed . 
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Introduction 


The  Conservation  Needs  Inventory  (CNI)  of  1967  is  an  important  data  base 
for  soil,  land  use  and  treatment  needs  information  which  is  consistent 
throughout  the  United  States.  In  the  past,  however,  CNI  information  has  been 
used  in  an  aggregated  form,  primarily  as  tables  of  land  use  and  soils  by 
county  or  region.  While  the  sampling  method  employed  by  the  CNI  is  only 
valid  for  individual  counties  or  larger  geographic  areas,  much  of  the  spatial 
information  contained  in  the  data  has  not  been  used. 

This  paper  explains  two  methodologies  which  make  fuller  use  of  the 
spatial  dimension  inherent  in  the  individual  point  samples  of  the  CNI  data 
base:  1)  computer  mapping  based  on  the  CNI  primary  sample  units  (PSUs);  and 
2)  computation  of  gross  soil  erosion  using  the  Universal  Soil  Loss  equation 
evaluated  at  CNI  sample  points.  Applications  to  the  Chowan-Pasquotank  River 
basin  of  Virginia  and  North  Carolina  are  presented.  Listings  of  the 
computer  programs  developed  are  included  in  appendices. 

Such  techniques  should  increase  the  usefulness  of  the  existing  CNI 
data  base  as  well  as  the  1977  update.  Computer  methods  used  in  conjunction 
with  this  data  base  can  result  in  more  thorough,  varied  analysis  of  river 
basin,  regional  and  other  planning  issues.  Such  methods  are  likely  to  be 
more  economical  as  well. 
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Planning  Limitations  of  Existing  Soil  Maps 


Mapped  soils  information  is  generally  available  to  planners  in 
two  forms:  1)  as  detailed  soil  maps  showing  the  location  of  each  soil 
on  a  relatively  large  scale  (small  area)  map;  or  2)  as  soil  association 
maps  showing  the  general  location  of  broad  groups  of  soils.  Detailed 
soil  maps  are  often  not  useful  for  planning  efforts  at  the  river  basin 
level  because  they  provide  an  overabundance  of  information.  Soil  associa¬ 
tion  maps,  on  the  other  hand,  obscure  important  distinctions  relevant  to 
large  area  planning  and  are  often  not  consistent  across  political 
boundari es . 

Soil  associations  are  logical  groupings  of  soils  with  similar 
parent  material  that  are  found  physically  to  occur  together.  Although 
associated  geographical ly ,  the  soils  in  an  association  may  differ  in 
their  suitability  for  agricultural  production.  For  example,  the 
Norfol k-Wagram-Rai ns-Bi bb  soil  association  in  Northampton  County,  NC, 
shown  in  Figure  1,  consists  of  soils  which  are  good  (60%),  good  to 
fair  (20%),  fair  (10%),  and  poor  (10%)  for  agricultural  use. 

A  second  difficulty  with  some  of  the  more  detailed  soil  association 
maps  is  that  they  were  originally  prepared  for  Soil  and  Water  Conservation 
Districts,  which  approximate  county  boundaries.  Different  districts 
group  soils  into  different  associations,  making  cross-county  comparisons 
in  larger  planning  areas  difficult.  The  problem  is  compounded  in 
multi-state  planning  areas,  as  shown  in  Figure  1,  because  not  only 
are  associations  composed  of  different  soils,  but  similar  soils  are  often 
given  different  names  in  different  states.  These  problems  can  be 
overcome  only  through  considerable  time  and  effort  on  the  part  of  soils 
experts . 
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Detailed  soil  maps  avoid  the  difficulties  encountered  in  using 
soil  associati ons because  the  productivity  characteri sties  of  individual 
soils  are  known  and  because  they  are  mapped  as  discrete  units.  An 
overaboundance  of  information  makes  detailed  soil  maps  less  useful  for 
planning  purposes.  For  example,  there  are  18  different  soils  mapped  in 
just  one  square  mile  of  Figure  2.  Only  about  55%  of  the  nation  has  been 
mapped,  so  many  areas  across  the  country  do  not  have  detailed  soil  maps. 
Even  when  available,  detailed  soil  maps  must  be  generalized  to  a  more 
useable  level  of  detail.  Also,  river  basin  planning  studies  produce  large 
amounts  of  information  which  cannot  be  effectively  overlaid  with  detailed 
soil  maps  without  rescaling.  This  paper  will  argue  that  computerized 
soils  data  currently  available  can  be  used  much  more  efficiently  than 
detailed  soil  maps  that  would  first  have  to  be  digitized  or  handled  through 
cartographic  centers  for  generalization  and  rescaling. 

The  Conservation  Needs  Inventory:  A  Primer 

In  order  to  more  fully  appreciate  the  kind  of  information  available 
from  the  detailed  sample  records  collected  during  the  1967  Conservation 
Needs  Inventory  (CNI),  a  brief  review  of  the  methodology  and  terminology 
used  is  presented. 

Stratified,  random  sample  areas  were  selected  by  the  Statistical 
Laboratory  at  Ames,  Iowa,  by  dividing  each  county  in  the  nation  into  strata 
with  48  sample  areas  in  each.  One  sample  area  was  selected  at  random  from 
each  stratum.  Sample  areas  varied  from  40  to  640  acres  each  and  were 
mapped  on  county  base  maps  supplied  by  the  State  SCS.  These  sample  areas 
are  called  Primary  Sample  Units  (PSU)  and  represent  approximately  two 
percent  of  the  area  of  each  county.  (Soils  Memorandum  21 (revision  2), 
November  10,  1965). 


4 


P-Smyra'; 

Church'', 


LUMffCR  «u.  •  j 

X: \  '■  •  ■rV'T'^wJ1 


Brn#V' 


Figure  2--Detailed  Soil  Map,  Northampton  County,  NC 
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The  PSUs  were  superimposed  on  soil  maps  of  the  county  and  a  sampling 
template  oriented  randomly  over  the  PSU  to  obtain  several  dozen  individual 
sample  points.  Soil  information  for  each  of  the  individual  sample  points 
was  read  directly  from  the  soil  map  and  land  use  and  treatment  information 
was  obtained  through  field  checks.  (Conservation  Needs  Memorandum  21, 
November  3 ,  1 965) . 

The  Statistical  Laboratory  expanded  the  PSU  data  to  inventory  acreages 
of  land  capability  classes  and  subclasses  obtained  during  the  same  year, 
1967.  Data  was  reviewed  and  adjusted  by  state  and  county  CNI  committees 
organized  by  the  State  SCS. 

Magnetic  tapes  were  prepared  for  each  state  containing  one  record  of 
24  fields  for  each  individual  sample  point.  Each  record  has  state,  county, 
soil  conservation  district,  watershed  and  PSU  identifying  codes;  soil, 
slope,  and  land  capability  codes;  land  use  and  conservation  treatment  need 
codes;  and  an  acreage  showing  the  area  of  the  county  represented  by  the 
point.  There  are  149,459  records  for  North  Carolina  and  97,980  for 
Virginia,  with  36,505  records  in  the  26  counties  of  the  Chowan  Basin. 

A  subsample  of  the  1967  CNI  was  resampled  in  1977.  Two  to  three 
points  for  selected  PSUs  were  sampled  obtaining  updated  land  use  and 
treatment  needs,  data  necessary  to  evaluate  the  Universal  Soil  Loss 
Equation  (USLE)  and  other  data. 

An  Alternative  Method  for  Mapping  Soils  Information 

Soils  data  for  broad  planning  purposes  can  be  mapped  relatively  easily 
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using  computer  based  data.  Soil  mapping  units  thus  derived  are  consistent 
across  county  and  state  boundaries,  like  detailed  soil  maps,  but  soils 
are  grouped  for  planning  purposes.  The  process  entails:  1)  grouping  soils 
into  units  meaningful  for  the  particular  study;  2)  sampling  the  location 
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of  these  units  in  space;  and  3)  mapping  the  sampled  locations  and  their 
associated  values.  A  description  of  the  three  parts  of  the  process 
fol 1 ows . 

Establishing  Meaningful  Soil  Groups  -  The  first  step  in  the  process 
is  to  group  individual  soils  into  fewer  units  which  will  be  easier  to 
manipulate  and  which  will  have  characteristics  meaningful  to  the  study. 

For  the  Chowan  River  Basin,  agricultural  productivity  was  of  primary  inter¬ 
est,  hence  soil  productivity  groups  (SPG's)  were  identified  and  created. 

The  soils  of  North  Carolina  were  arranged  into  32  groups  based  on  crop 
yield  characteri sties ,  management  practices,  and  land  treatment  needs 
for  the  Tar-Neuse  River  Basin  study  using  CNI  data.  (Parks,  1975). 

Because  Virginia  soils  names  differ  slightly  from  those  used  for  the  same 
soils  in  North  Carolina,  the  help  of  Virginia  soil  scientist,  Richard 
Googins,  was  enlisted  to  extend  the  North  Carolina  SPG's  to  Virginia 
soils.  Thus,  a  uniform,  consistent  classification  of  soils  based  on 
productivity  and  conservation  needs  was  established  transcending  county 
and  state  boundaries. 

Locating  Sample  Units  -  In  order  to  map  SPG's  it  is  necessary  to 
determine  their  location  in  space  at  a  large  number  of  points.  For  this 
purpose  the  CNI  sample  record  provided  an  adequate  sample  of  the  location 
of  SPGs . 

Commensurate  with  the  level  of  detail  necessary  in  the  data  and  in 
order  to  reduce  the  required  effort,  the  36,505  CNI  records  within  the 
Chowan  Basin  were  aggregated  to  the  PSU  level  by  selecting  the  major 
SPG,  land  use  and  treatment  need,  on  the  basis  of  the  acreage  represented 
on  the  individual  records  making  up  the  PSU.  The  process,  diagrammed  in 
Figure  3,  resulted  in  1,408  aggregated  records,  OTie  record  for  each  PSU, 
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Figure  3 
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for  a  96%  reduction  in  data  from  the  original  36,505  sample  point  records. 

Since  CNI  records  contain  no  geographical  coordinates,  microfilm 
copies  of  the  county  base  maps,  with  PSU's  indicated,  were  obtained  and 
large  paper  copies  produced  on  which  latitude  and  longitude  lines  were 
superimposed.  The  coordinates  of  each  PSU  were  determined  to  the  nearest 
half  minute  and  were  appended  to  the  aggregated  CNI  record  for  each  PSU. 

The  resulting  data  base  consisted  of  1,408  records,  one  for  each  PSU,  with 
the  dominant  SPG,  1967  land  use  and  1967  treatment  need,  and  geographic 
coordinates.  Including  geographic  coordinates  among  the  data  contained 
in  each  sample  point  record  would  greatly  improve  the  CNI  and  any  such 
future  data  gathering  efforts. 

Mapping  Sample  Units  -  In  order  to  graphically  display  the  data,  coor¬ 
dinates  on  the  aggregated  PSU  records  were  combined  with  county  boundaries 
contained  on  the  Census  Bureau's  Dual  Independent  Map  Encoding  County 
File  (DIMECO)  as  shown  in  Figure  4. 

First,  the  county  boundaries  of  interest  in  the  Chowan  Basin  were 
selected  from  the  DIMECO  file,  which  contains  boundaries  for  all  states 
and  counties.  Next,  the  segment  file  structure  of  DIMECO  was  converted 
into  a  different  file  structure  in  which  "chains"  of  segments  making 
up  entire  polygons  were  strung  together  from  the  individual  boundary 
segments.  This  is  an  important  step  in  that  most  mapping  program  packages, 
such  as  SYMAP  and  CALFORM,  use  a  form  of  this  chain  file  structure  in  their 
input  data.  Conversion  of  the  Basin  boundary  file  was  accomplished  using 
POLYVRT,  a  program  package  produced  by  the  Harvard  Laboratory  for  Computer 
Graphics  and  Spatial  Analysis.  (POLYVRT  User's  Manual,  November,  1974) 

Two  kinds  of  graphical  display  were  produced  from  the  data  contained 
in  the  aggregated  PSU  records  and  in  the  Basin  boundary  file:  1)  plotter 
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Figure  4 
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maps  showing  the  boundary  of  counties  in  the  Basin  and  the  location  and 
identity  of  each  PSU,  as  in  Figure  5;  and  2)  printer  maps  showing  shaded 
zones  based  on  the  PSU's  as  in  Figure  9.  The  programs  which  accomplish 
this  mapping  are  briefly  described  below  and  are  discussed  in  more 
detail  in  Appendices  A-l  and  A-2. 

Plotter  mapping  is  the  simpler  of  the  two  forms,  and  is  accomplished 
by  a  program  called  CH0WMAP1  which  was  written  for  this  project.  The 
program  scales  the  geographic  coordinates  of  boundary  points  and  data 
points,  converts  them  to  a  Cartesian  coordinate  system  used  by  the  plotter, 
limits  the  plot  to  a  "window"  specified  by  control  cards  through  a 
process  called  "clipping,"  and  controls  the  movements  of  the  plotter 
pen  to  draw  the  county  boundaries  and  indicate  the  location  and  identification 
of  each  data  point.  Data  points  consisted  of  the  PSU's,  which  could  be 
identified  by  the  associated  SPG,  land  use,  or  treatment  need,  or  by  any 
other  identification  field  in  the  aggregated  PSU  record. 

Printer  mapping  is  more  complicated.  As  shown  in  Figure  6,  the 
process  begins  with  the  initial  map  matrix  on  which  land  boundaries  and 
the  SPG  identifications  for  each  PSU  are  located.  For  each  empty  cell 
in  the  map  matrix  which  is  inside  the  land  boundaries.,  an  SPG  is  assigned 
based  on  the  SPG  of  the  nearest  PSU.  When  all  empty  cells  within  the 
boundary  are  assigned,  the  map  matrix  is  complete.  The  computer 
program  actually  works  with  a  numerical  representation  of  the  map  matrix 
called  a  scandeck  (See  Appendix  A-2).  Using  the  computer  mapping  package 
SCANGEN  (Center  for  Census  Use  Studies,  May,  1975)  any  vector  of  data 
associated  with  SPG's  can  be  mapped  into  the  zones  defined  by  the  scandeck. 
Thus,  zones  are  created  only  once  and  stored  for  use  in  mapping  various 
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Figure  5 


Plotter  map  of  PSU's  identified  by  Soil  Productivity  Group  (SPG) 
Albemarle  Sound,  North  Carolina 


Figure  6--  Creation  of  Scandeck  and  Printer  Maps 
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data  associated  with  them.  This  is  a  distinct  advantage  over  similar 
mapping  packages  which  must  create  zones  each  time  a  map  is  desired. 

Application  to  the  Chowan  Basin 

Tv/o  applications  of  this  methodology  to  river  basin  planning  in  the 
Chowan-Pasquotank  Basin  will  be  shown  to  illustrate  the  utility  of  the 
approach.  First,  the  amount  and  location  of  land  use  change  around 
Albemarle  Sound  was  investigated.  Second,  the  amount  and  location  of 
gross  soil  erosion  from  agricultural  activities  in  the  Basin  was  calcula¬ 
ted  and  mapped. 

Amount  and  Location  of  Land  Use  Change  -  Establishing  which  Soil 
Productivity  Groups  are  experiencing  rapid  change  to  agricultural  uses, 
and  which  land  uses  are  most  affected,  indicates  the  potential  for  future 
change  through  land  clearing. 

The  aggregated  PSU  records  contain  the  dominant  land  use,  in  terms  of 
area,  of  the  PSU  from  the  1967  Conservation  Needs  Inventory  (CNI).  Rapid 
clearing  and  draining  of  forested  wetlands  made  CNI  land  use  data  obsolete. 

In  order  to  update  the  land  use  information  for  wetland  soils  two  kinds  of 
more  recent  land  use  data  were  obtained.  The  United  States  Geological  Survey 
publishes  rectified  7.5  minute  orthophoto  quadrangle  maps  at  a  scale  of 
1:24,000  which  display  adequate  shading  and  texture  to  distinguish  between 
major  land  uses  wuch  as  cropland,  pasture,  forest,  and  urban  areas.  In 
addition,  for  more  limited  areas,  USGS  has  prepared  Land  Use  and  Land  Cover 
maps  from  airphoto  and  satellite  data  at  a  scale  of  1:250,000.  (Anderson, 
et  al.,  1976)  These  maps  distinguish  between  37  different  land  uses  and 
cover  types,  shown  in  Table  1,  and  are  thus  more  detailed  than  the 
orthophoto  maps.  Land  use  and  land  cover  maps  showing  1974  land  use 
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were  available  from  USGS  for  Rocky  Mount  and  Manteo  quadrangles. 
Orthophoto  maps  showing  1974  land  use  were  obtained  for  most  of  the 
rest  of  the  region  surrounding  Albemarle  Sound  from  the  North  Carolina 
Department  of  Natural  and  Economic  Resources,  as  shown  in  Figure  7. 

Land  use  change  between  1967  and  1974  was  assessed  using  the  ortho¬ 
photo  quad  maps  in  the  Northern  Albemarle  Sound  area  by  first  super¬ 
imposing  a  ^-minute  grid  over  each  quadrangle  and  locating  each  PSU 
by  means  of  the  coordinates  on  each  aggregated  PSU  record.  The  1974 
land  use  was  recorded  as  one  of  four  possibilities:  1)  cropland  and 
pasture;  2)  forest  land;  3)  urban  land;  or  4)  water.  Each  PSU  record 
thus  processed  contained  a  land  use  observation  for  two  points  in 
time.  Since  the  land  use  and  land  cover  maps  were  available  as  trans¬ 
lucent  overlays,  a  plot  map  of  the  county  outlines  and  aggregated  PSU 
data  points  for  the  Southern  Albemarle  Sound  area  was  created  at  the 
same  scale  as  the  land  use  overlay  using  CHOWMAP  1  and  the  1974  land 
use  read  directly  from  the  combined  map. 

The  1974  land  use  data  from  both  sources  was  appended  to  the 
aggregated  PSU  records  so  that  1967  land  use,  1974  land  use,  or  the 
change  in  land  use  from  1967  to  1974  could  be  determined  and  mapped  using 
CHOWMAP  2.  This  data  could  also  be  displayed  in  tabular  form,  as  in 
Table  2,  by  direct  processing  of  the  aggregated  PSU  records. 

Amount  and  Location  of  Gross  Soil  Erosion  -  Data  from  individual 
CNI  sample  records  was  used  in  conjunction  with  the  Universal  Soil  Loss 
Equation  (USLE)  to  project  1990  gross  soil  erosion  by  Soil  Productivity 
Group  in  the  Chowan-Pasquotank  Basin.  Erosion  was  projected  under 
different  scenarios  as  to  future  production  levels  in  the  Basin,  and 
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Table  1--Land  use  and  land  cover  classifications 
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LAND  USE  AND  LAND  COVER  MAPS  The 

LEVEL  1 

1  Urban  or  Built-up  Land 


2  Agricultural  Land 


3  Rangeland 


4  Forest  Land 


5  Water 


6  Wetland 

7  Barren  Land 


8  Tundra 


9  Perennial  Snow  or  Ice 


mans  show  the  following  categories. 
LEVEL  II 
Resident  1  a l 

Commercial  ami  Services 
I nd us  tria  1 

Transportation ,  Communications , 
and  Utilities 

Industrial  and  Commercial  Complexes 
Mixed  Urban  or  Built-up  Land 
Other  Urban  or  Built-up  Land 

Cropland  and  Pasture 

Orchards,  Groves,  Vineyards,  Nurserie 
and  Ornamental  Horticultural  Areas 
Confined  Feeding  Operations 
Other  Agricultural  Land 

Herbaceous  Rangeland 
Shrub  and  Brush  Rangeland 
Mixed  Rangeland 

Deciduous  Forest  Land 
Evergreen  Forest  Land 
Mixed  Forest  Land 

Streams  and  Canals 
Lakes 

Reservoirs 

Bays  and  Estuaries 

Forested  Wetland 
Non fores  ted  Wetland 

Dry  Salt  Flats 
Beaches 

Sandy  Areas  other  than  Beaches 
Bare  Exposed  Rock 

Strip  Mines,  Quarries,  and  Gravel  Pit 
Transitional  Areas 
Mixed  Barren  Land 

Shrub  and  Brush  Tundra 
Herbaceous  Tundra 
Bare  Ground  Tundra 
Wet  Tundra 
Mixed  Tundra 

Perennial  Snowfields 
Glaciers 


se 

I  l 

I  2 

L  3 

14 

1  5 

16 

17 

21 

22 

23 

24 

31 

32 

33 

41 

42 

43 

51 

52 

53 

54 

61 

62 

71 

73 

73 

74 

75 

76 

77 

81 

82 

83 

84 

85 

91 

92 
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Table  2--  Soil  Productivity  Groups  with  Land  Use  Changes, 

Northern  Albemarle  Sound- 
1967-1974 


PG 

1967 

Cropland 
&  Pasture 

New 

+  Cropland  & 

Pasture  Uses 

New 

Forest 

Uses 

1974 

Cropland 
&  Pasture 

1967-1974 

Net 

Change 

2 

22,628 

0 

4,056 

18,572 

-4,056 

3 

60,303 

23,340 

18,423 

65,220 

4,917 

5 

38,518 

28,304 

11,248 

55,574 

17,056 

5 

133,476 

39,720 

22,521 

150,675 

17,199 

1 

23,891 

24,833 

6,724 

42,000 

18,109 

i 

0 

1,292 

0 

1,292 

1,292 

107 

11,005 

0 

11,112 

11,005 

)TAL 

278,923 

128,494 

62,972 

344,445 

65,522 

1/ 

Camden,  Chowan, 

Currituck,  Gates, 

Pasquotank,  and 

Perquimans  Counties. 

economic  effects  of  implementing  certain  soil  conservation  plans  were 
compared.  These  effects  include  changes  in  the  location  of  production, 
based  on  a  linear  programming  model  which  was  used  to  anticipate  farmers 
efforts  to  find  the  least  costly  way  of  meeting  environmental  objectives 

The  Universal  Soil  Loss  Equation  mul tipi icatively  combines  a  number 
of  factors  relevant  to  the  detachment  and  removal  of  soil  particles  on 
agricultural  land.  (Wischmeier  and  Smith,  1965)  The  equation  is: 

A  =  RxKxLxSxCxP 
where  the  factors  are  defined  as  follows. 

The  computed  soil  loss  per  acre.  A,  is  an  estimate  of  the  gross 
amount  of  soil  which  is  detached  from  the  land  surface  over  the  course 
of  a  crop  year.  No  account  is  made  for  deposition  of  sediment  on  other 
parts  of  the  same  field,  so  it  does  not  measure  the  escape  of  sediment 
from  the  field. 

The  rainfall  factor,  R,  is  derived  from  a  map  of  iso-erodent  values 

measured  empirically.  The  R  factor,  like  many  of  the  USLE  factors,  is 
tabulated  in  the  North  Carolina  Soil  Conservation  Service  Technical 
Guide.  (NCSCS ,  October,  1976)  The  computer  program  which  calculates 
soil  loss,  diagrammed  in  Figure  8,  selects  the  R  factor  value  based  on 
the  county  code  from  each  CNI  sample  record  from  a  list  of  R  factors  tab 
ulated  by  county,  taken  from  the  SCS  Guide.  The  soil  erodability 
factor,  K,  is  tabulated  by  soil  type  in  the  SCS  Guide.  Examination  of 
soils  data  for  SPG's  in  the  Chowan  Basin  showed  that  for  each  SPG 
fewer  than  four  soils  always  accounted  for  more  than  two-thirds  of 
the  SPG  area.  The  EROSION  (Appendix  A-3)  program  examines  the  soil  code 
on  each  individual  CNI  sample  record,  determines  if  it  is  a  dominant 
soil,  and  assigns  the  correct  K  factor  value. 
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Figure  8— Soil  loss  computation  flow  chart 


Modify  L  value 
and  recompute 
RKLS 


Select  F  values 
and  compute 
RKLSP 


Print 


Go)  to  1 
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The  slope  length  factor,  L,  is  perhaps  the  most  difficult  factor 
to  evaluate  given  present  data.  A  survey  conducted  in  1972  by  the  Soil 
Conservation  Service  determined  the  dominant  soil,  slope  length,  slope, 

K  value  and  T  value  for  land  resource  areas  in  selected  states.  (McElroy, 
et  al..  May  1976,  and  correspondence  with  Roy  M.  Gray,  SCS  dated  4  April 
1977)  Using  this  data  a  slope  length  was  related  to  each  dominant  soil. 
Thus,  a  slope  length  is  automatically  assigned  when  the  K  value  is 
determi ned . 

The  slope  gradient  factor,  S,  was  taken  directly  from  the  slope 
field  of  the  individual  CNI  sample  record.  The  slope  gradient  and  slope 
length  values  are  used  in  the  equation  for  the  LS  factor  contained  in 
SCS  (September,  1977,  shown  in  Figure  8. 

The  product  RKLS  is  computed  from  the  factors  listed  above  and 
constitutes  the  basic  soil  erosion  value  associated  with  a  particular 
sample  point  in  a  particular  SPG,  listed  in  Table  3.  The  RKLS  values 
for  each  sample  point  in  an  SPG  are  averaged,  using  the  acreages  repre¬ 
sented  by  each  sample  point  as  weights  to  obtain  an  area-weighted  average 
bare  ground  soil  loss  for  the  SPG.  This  value  is  later  multiplied  by 
the  cropping  management  factor,  C,  for  the  appropriate  crops  and 
rotations  and  by  the  erosion  control  practice  factor,  P,  which  are 
tabulated  in  the  SCS  Guide.  Thus,  an  estimate  of  gross  soil  erosion  for 
the  particular  soils,  slopes,  cropping  practices,  and  conservation 
practices  is  calculated.  Plans  for  conserving  soil  and  preventing  nonpoint 
pollution  primarily  analyze  the  economic  and  locational  impacts  of 
practices  with  low  P  and  C  values  on  crops  and  soils  with  erosion  potential. 
The  results  of  one  such  plan  are  shown  in  Figure  9  for  the  Chowan  Basin. 
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Table--  3 

Gross 

soil  erosion,  Chowan  Basin 

1/ 

SPG. 

No 

Practice 

Contour 

Farming 

Field 
SC  A 

2/ 

Field 

SC  B  3/ 

Contour 

SC  A  4/ 

Contour 
SC  R 

10. 

61.2 

30.6 

21.4  24.5 

15.3 

21.4 

11. 

95.6 

57.4 

38.2 

43.0 

28.7 

3S .  2 

12. 

75.4 

37.7 

26.4 

30.2 

18.9 

26.4 

13. 

29.2 

17.5 

10.2 

13.1 

8.7 

11." 

14. 

98.1 

49.1 

34.3 

39.2 

24.5 

34 . 3 

15. 

18.5 

11.1 

6.5 

8. 3 

5.6 

7.4 

16. 

138.9 

83.4 

55.6 

62.5 

41.7 

55.6 

17. 

131.1 

104.9 

65.6 

85.2 

52.5 

65.6 

18. 

41.2 

22.1 

14.4 

17.2 

11.0 

15.1 

19. 

106.4 

53.7 

37.5 

42.8 

26.9 

37.5 

20. 

127.2 

63.6 

44.5 

50.9 

31.8 

44.5 

21. 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

22. 

31.1 

17.5 

10.9 

13.4 

8.7 

11. S 

23. 

41.1 

24.7 

14.4 

18.5 

12.4 

16.5 

24. 

29.1 

17.5 

10.2 

13.1 

8.7 

11.7 

25. 

20.5 

12.3 

7.2 

9.2 

6.1 

8 . 2 

26. 

37.3 

22.4 

13.0 

16.7 

11.2 

14.9 

27. 

47.0 

28.2 

16.4 

21.1 

14.  1 

18.8 

28. 

20.8 

12.4 

7.2 

9.4 

6.3 

8.3 

29. 

167.0 

94.0 

63.7 

75.5 

47.0 

65 . 7 

30. 

42.6 

22.1 

14.9 

17.4 

11.0 

15.3 

31. 

1.1 

0.7 

0.4 

0.5 

0.3 

0.5 

\J  Soil  losses  shown 
to  bare-ground  or 

do  not  include  the  effects 
fallow  soil  lost. 

of  the  crop 

practice 

factor,  C, 

and 

are  thus 

equivalent 

2/  Field 

stripcropping  with  meadow  strips. 

V  Field 

stripcropping  with  small 

grain  winter  < 

:over . 

-  c°ntour  stripcropping  with  meadow  strips. 

j  Contour  stripcropping  with  small  grain  winter  cover. 
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Figure  9 --Location  of  cropland  with 

projected  for 


various 

1990 


erosion  rates 
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This  method  for  obtaining  gross  soil  erosion  can  be  favorably  com¬ 
pared  with  similar  methods  which  also  use  the  Universal  Soil  Loss  Equation.) 
Factors  are  evaluated  for  each  sample  point  and  the  soil  losses  averaged, 
rather  than  averaging  the  factor  values.  Figure  10  compares  the  method 
used  in  this  study  with  two  alternative  methods  for  computing  soil  loss. 

The  first  method  involves  averaging  the  factor  values  of  the  USLE  over 
the  entire  Basin,  a  process  which  naively  overlooks  huge  differences  in 
conditions  over  an  area  this  large.  Method  2  ameliorates  these  differences 
somewhat  by  averaging  the  USLE  factors  for  sub-basins  judged  to  be  more 
homogeneous.  The  third  method,  used  in  this  study,  computes  the  soil 
loss  for  many  points  within  the  Basin  and  averages  the  soil  loss  for  the 
soils  making  up  each  SPG  rather  than  averaging  the  factor  values.  This 
method  more  completely  accounts  for  differences  in  rainfall,  soil,  and 
slope  over  the  Basin  and  is  likely  to  be  more  accurate.  For  making  current 
soil  loss  estimates  or  projections,  this  method  should  be  even  more 
acceptable  when  the  results  of  the  1977  Soil  Conservation  Service  Erosion 
Inventory  are  available.  (SCS,  June,  1977)  This  inventory  resurveys 
a  subset  of  the  original  sample  points  consisting  of  two  to  three  sample 
points  within  specified  PSU's.  Land  use  and  treatment  information  will 
be  updated  and  information  on  specific  types  of  crops  grown,  conserva¬ 
tion  practices  used,  and  USLE  factors,  including  slope  length  factors 
which  were  missing  from  the  1967  inventory,  will  be  added. 
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Figure  10--  Three  Methods  for  Computing  Soil  Loss,  Chowan  Basin 
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Concl usion 


The  methods  presented  in  this  paper,  in  conjunction  with  the  CNI 
sample  data,  offer  a  spatial  dimension  to  agricultural  resource  analysis 
which  has  largely  been  absent  in  the  past.  The  ability  to  determine  and 
spatially  display  resource  attributes  opens  up  a  wide  range  of  possible 
applications  to  river  basin  planning  and  to  other  studies. 

For  river  basin  planning,  these  methods  aid  both  the  creation  of 
plans  and  their  evaluation.  Identification  of  changing  agricultural  areas, 
in  terms  of  new  clearing  or  in  terms  of  urban  encroachment,  can  be  made 
by  comparing  1967  land  use  from  the  CNI  to  later  land  use  data  from  a 
variety  of  sources.  Analysis  of  soils  on  which  change  occurs  can  further 
focus  planning  efforts  on  those  areas  important  to  a  basin's  future. 

Impacts  of  alternative  plans  can  be  displayed  graphically  for  easy  compari¬ 
son.  The  mapping  programs  presented  here  are  designed  to  display  any  of 
27  outputs  from  an  LP  model,  located  by  means  of  SPG.  These  outputs  include 
economic  information  such  as  cost  of  agricultural  production,  practice 
information  such  as  amount  of  terracing,  production  information  on  each 
crop  grown  and  environmental  information  such  as  amounts  of  fertilizer 
applied.  The  output  of  any  other  analysis  by  soil  productivity  groups  can 
also  be  displayed.  Because  the  method  used  stores  the  zone  information 
created  by  the  nearest-neighbor  technique,  each  subsequent  map  is  inexpen¬ 
sive  to  create.  Changes  in  shading  or  interval  values  can  be  accommodated 
easily. 
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Possible  applications  exist  for  other  studies  as  well.  Water 
quality  studies  can  benefit  from  a  more  precise  computation  of  soil 
erosion,  which  is  of  importance  in  nitrogen  and  phosphorus  movement  as 
well  as  in  terms  of  sediments.  Analysis  of  farm  production  policies 
can  benefit  from  spatial  displays  of  their  potential  impact  on  the 
location  of  production.  The  cost  and  ease  of  application  of  these  methods 
open  up  the  possibility  of  conducting  fairly  detailed  planning  work  at 
the  state  or  regional  level.  This  could  benefit  water  quality,  conserva¬ 
tion  arid  pesticide  studies  which  have  relatively  short  turn-around  and 
staffing.  Other  analyses  will  likely  be  possible  in  the  future,  to  the 
extent  that  update  studies  continue  to  build  on  the  existing  CNI  data 
base,  adding  new  and  different  data  items.  The  potential  for  the  CNI 
to  become  the  basis  for  a  national  land  and  water  resource  data  base 
certainly  exists  and  is  enhanced  by  development  of  these  and  other 
methods  to  utilize  such  a  data  base. 
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Appendix  A-l 


CH0WMAP1 

A  computer  mapping  program  that  produces  county  outline  maps 
with  PSU's  identified  by  SPG  on  the  CALCOMP  plotter. 

Inputs- 

Unit  8  U.P006066. CHPWAN.POLYVRT.CHAI NS- -chain  file  of  county 

boundary  segment  chains. 

Uni t  9  U . P006066 . CHOWANNC . CNI COORD. PSUSORT 
U. P006066 . CHOWANVA. CNICOORD. PSUSORT 
coordinates  and  identifying  data  for  PSU's  in  North  Carolina  and 
Virginia  portions  of  the  Chowan  Basin. 

Outputs- 

Description  of  area  mapped,  window  of  map  and  points  plotted. 
Plotter  map  (see  Figure  5). 

Run  Time- 

Two  quadrangle  maps 
Rocky  Mount,  NC 

172  vertices  in  outline 
165  data  points  (PSUs) 

Northern  Albemarle  Sound,  NC 

465  vertices  in  outline 
830  data  points  (PSUs) 

8.0  seconds  of  IBM  370/168  CPU 
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Appendix  A-2 


CH0WMAP2 

A  collection  of  computer  programs  that  produces  chloropleth  maps 
of  zones  based  around  PSUs  on  the  line  printer.  The  programs  are: 
MATMAP--produces  a  numerical  representation  of  the  map  matrix  called 
a  scandeck  which  can  be  stored  for  reuse. 

Inputs- 

Unit  5  Control  cards  specifying  the  kind  of  map  to  be  produced, 
the  geographic  coordinates,  the  size  of  the  viewpoint,  the  character¬ 
istics  of  the  printer,  a  map  title  and  the  input  file  format. 

Unit  8  Description  of  the  input  file  including  number  of  chains, 
number  of  polygons,  number  of  vertices,  and  the  number  of  nodes.  The 
input  file  itself  is  also  read  on  this  unit  and  consists  of  chains 
of  boundary  segments. 

Unit  9  Input  file  of  data  points  which  will  form  the  zone 
centroids. 

Outputs- 

Unit  FTOUT  Output  file  for  scandeck.  An  example  of  a  scandeck 
is  included  in  this  appendix. 

REFORMAT  --  Reformats  data  to  be  mapped  and  selects  PSU  coordinates 
for  use  in  SCANGEN. 

Inputs- 

Unit  1  Data  to  be  mapped,  by  SPG 
Unit  5  Coordinates  of  map  frame 
Unit  8  Coordinates  of  PSUs 
Outputs- 

Unit  9  Reformated  file  containing  PSU  coordinates,  SPG  associated 
with  the  PSU  and  data  value  to  be  mapped  associated  with  that  PSU 
through  the  SPG. 

SCANGEN  (Center  for  Census  Use  Studies,  May,  1976 )--Actual ly  produces 
the  map. 

Inputs- 

Control  cards  specifying  input  file  and  output  map  characteristics. 
Outputs- 

A  description  of  the  map  produced  including  a  statistical  summary 
of  the  mapped  intervals  and  the  printer  map  (see  Figure  9). 
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Example  of  scandeck  output--each  line  represents  one  row  of  the  map 
matrix,  indicated  by  the  first  number.  The  number  of  segments  in 
the  line  (corresponding  to  the  number  of  zones,  including  background 
areas,  moving  from  left  to  right  across  the  map  matrix)  is  indicated 
by  the  second  number.  Pairs  of  numbers  following  this  give  the  zone 
identification  (internal  to  the  program)  and  the  ending  print  position 
for  the  zone. 

Lines  1  through  13  below  are  all  background  (water)  lines  con¬ 
sisting  of  1  zone  identified  as  958  which  ends  in  print  position  132. 
Line  14  is  the  first  line  in  which  some  mapped  area  of  interest  shows 
up.  It  has  2  zones,  the  first  identified  as  zone  4  which  ends  in 
print  position  4,  and  the  second  which  is  background. 


1  1 9581 32 

2  1958132 

3  1  9  381  32 
A  IV 1381  32 

5  1958132 

6  1 7  3  o 1 32 

7  IV 381 32 

8  1 V  38 1  32 

9  1958132 

10  1938132 

11  1938132 
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Appendix  A-3 


EROSION 

A  computer  program  which  calculatres  gross  bare-ground  soil 
loss  using  the  universal  soil  loss  equation  with  factors  evaluated 
for  each  sample  point  in  the  Conservation  Needs  Inventory. 

Inputs- 

Unit  10  *. SORT. SORT. S0RT0UT--i ndi vi dual  CNI  sample  records 
sorted  by  soil  productivity  group. 

Unit  11  U. P006066. RH. EROSION. DATA- -Factor  values  of  the  USLE. 

R  factor  by  county 

P  factor  by  slope  interval  and  practice 
K  and  L  factors  by  dominant  soil 

Output- 

Table  of  normal  and  terraced  gross  bare-ground  soil  loss  estimates 
by  SPG  and  conservation  practice,  acreage  over  which  computed  and 
acreage-weighted  average  slope. 
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